Pair correlation function of short-ranged square-well fluids 
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We have performed extensive Monte Carlo simulations in the canonical (NVT) ensemble of the 
pair correlation function for square-well fluids with well widths A— 1 ranging from 0.1 to 1.0, in units 
of the diameter a of the particles. For each one of these widths, several densities p and temperatures 
T in the ranges 0.1 < pa 3 < 0.8 and T C (A) < T < 3T C (A), where T C (A) is the critical temperature, 
have been considered. The simulation data are used to examine the performance of two analytical 
theories in predicting the structure of these fluids: the perturbation theory proposed by Tang and 
Lu [Y. Tang and B. C.-Y. Lu, J. Chem. Phys. 100, 3079, 6665 (1994)] and the non-perturbative 
model proposed by two of us [S. B. Yuste and A. Santos, J. Chem. Phys. 101, 2355 (1994)]. It is 
observed that both theories complement each other, as the latter theory works well for short ranges 
and/or moderate densities, while the former theory does for long ranges and high densities. 
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I. INTRODUCTION 

The study of the thermodynamic and structural prop- 
erties of square-well (SW) fluids has been the subject 
of interest for many years because of their simplicity 
and their resemblance with real fluids with spherically 
symmetrical potentials, among other reasons. There- 
fore, at present there are available a considerable num- 
ber of theories for this kind of fluids. Among them, 
particularly simple and fruitful are perturbation theo- 
ries for the thermodynamic propertiesiSiiMiSiSiL&iLiil If 
one is interested in structural properties, one can resort 
to integral equations theories based on the Ornstein- 
Zcrnikc equation, for which we have a number of pos- 
sible choices, 11 ! 12 ! 13 ! 14 . 15 ! 16 ! 17 ! 18 ! 19 ! 20 ! 21 ! 22 ! 23 ! 24 The latter 
group of theories have in general the drawback of being 
non-analytical, so one has to deal with them by numer- 
ical methods. However, in some cases it has been pos- 
sible to obtain analytical expressions for the structural 



properties 



25,26.27,28.29,30 



inspired, at least indirectly, in 



data are used to test the performance of two analytical 
theories, one perturbativ o 27 ! 28 and the other one non- 
perturbativei 2 ^ As we will see, both theories complement 
each other: the perturbative theory is generally prefer- 
able for long ranges, while the non-perturbative theory 
is better for short ranges. 

The plan of the paper is as follows. The two theories 
are introduced in the next Section, some details being rel- 
egated to AppendiceslAlandlBl The Monte Carlo method 
we have employed is succinctly described in Section ITTT1 
The theoretical results are compared with the simulation 
data and discussed in Section llVI The main conclusions 
of the paper are summarized in Section IV! 



II. ANALYTICAL THEORIES FOR THE PAIR 
CORRELATION FUNCTION OF SQUARE- WELL 
FLUIDS 

For fluids with a square-well potential of the form 



integral equation theories. 

In parallel with the theoretical developments, 
a lot of research has been devoted to obtain- 
ing the thermodynamic and structural prop- 
erties of SW fluids by means of computer 
.:iimdati:iu:. R - in - 14 - 31 - 32 - 33 -^ 

Most of that research has focused on SW fluid with 
intermediate ranges of the potential, because they more 
closely mimic real simple fluids, whereas relatively little 
attention has been paid to SW fluids with short ranges. 
On the other hand, recently there has been a renewal 
in the interest in short-ranged SW fluids as models of 
colloidal suspensions^ 3 ^ 9 ^ * 5 ^^ and phase separation of 
protein solutions^ 

In the present paper we have carried out Monte Carlo 
simulations of the pair correlation function or radial dis- 
tribution function (r.d.f.) g(r) of SW fluids with short, 
intermediate, and long ranges for temperatures above the 
critical ones and for a wide range of densities. These 



. oo if r < a, 
u{r) = I —e if a < r < Act, 
if r > Act, 



(1) 



where A is the potential range in units of the particle di- 
ameter a and e is the potential depth, several approaches 
have been devised to derive analytical expressions for the 
structural properties. In this paper we will focus on two 
theories, both having in common that analytical expres- 
sions for the r.d.f. in Laplace space are provided. 

The first one of those theories is due to Tang and Lu 
(TL)j2L2£ who combined perturbation theory with the 
mean spherical approximation (MSA) to derive an an- 
alytical expression for the first-order r.d.f. gi(r) in the 
expansion in power series of the inverse of the reduced 
temperature T* — kT/e. Taking for the zeroth-order 
term go(r) the Percus-Yevick (PY) solution; 5 ^ 5 , the re- 
sulting truncated series for the r.d.f. of the SW fluid is 
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g(x)=g (x)+gi(x) — , (2) 

where x — r/a. The expressions for the Laplace trans- 
forms of xgo (x) and xg\ (x) are given in Appendix ^ 
The TL theory is expected to be accurate for moderate 
to large potential widths since the scries in powers of 
1/T* converges slowly for short-ranged SW potentials.- 6 

On the opposite situation, that is, for SW potential 
with ranges A close to 1, a procedure has been proposed 57 
to determine the structure of a SW fluid from that of 
an equivalent fluid of sticky hard spheres, using for the 
latter Baxter's analytical solution of the PY equation. 58 
To this end, the parameters of the equivalent fluid are 
determined from the condition that the second virial co- 
efficients of the two fluids must be equal. This approx- 
imation provides good results for the structure factor of 
SW fluids with A < 1.2, at least for moderate to low den- 
sities, but is not appropriate for obtaining the r.d.f., so 
that it will not be considered in this paper. 

As a second theory, we will consider the one devel- 
oped by Yuste and Santos? 2 ^ which provides an alter- 
native analytical expression for short-ranged SW fluids 
and reduces to Baxter's solution in the sticky hard-sphere 
limit. We will refer to this theory as the Yuste-Santos 
(YS) model and it will be presented next with some de- 
tail. 

The starting point in the YS model is the expression 
of the Laplace transform G(t) of xg{x) in the form 

G It) = t = y (-1277)"- 1 1 [F (t)] n e~ nt , 

w 1 + 12-nF (t) e~* ^ V " 1 K n 

' v / n=l 

(3) 

where rj — ^pcr 3 is the packing fraction, p being the 
number density, and F{t) is an auxiliary function given 

1 l + A + K.t-jA + K^e-^-^ 
{ ' 127 ? l + S 1 t + S 2 t 2 + S 3 t 3 ' 1 ' 

The coefficients Ki, K2, Si, S2, and S3 are determined 
from consistency conditions as functions of 77, T* , A, and 
A (see Appendix |B|) . To close the model, the parameter 
A is further fixed at its zero density value A = e 1 / 7 — 1 
for the sake of simplicityi22i22i Expression reduces to 
the exact solutions of the PY equation in the limit of 
hard spheres (A — * 1 or T* — > 00) j 5 ^ 5 . as well as in 
the limit of sticky hard spheres (A — > 1 and T* — > 
with T* ~ — 1/ ln(A— 1)),— Therefore, the approximation 
can be considered as an extension to finite widths 
of Baxter's solution of the PY equation for sticky hard 
spheres. 

The inverse Laplace transform of Eq. © allows us to 
obtain the r.d.f. in the form 

00 

g (x) = x- 1 (-Mr))"- 1 f n (x-n)Q(x- n), (5) 

71=1 



where the function /„ (x) is the inverse Laplace transform 
of t[F(t)] n and 0(x) is Heaviside's step function. Note 
that, to determine the r.d.f. for x < n+ 1 only the first n 
terms in the summation JSJ are needed. In the analysis 
of Section Hvl we will consider reduced distances x < 3, 
so that only the functions fi and f% will be needed. They 
are given in Appendix IBI 

III. MONTE CARLO SIMULATIONS 

We have performed NVT Monte Carlo (MC) simula- 
tions of the r.d.f. of SW fluids with ranges A = 1.05 and 
A = 1.1-2.0 (with a step AA = 0.1) for (reduced) number 
densities p* = pa z = 0.1-0.8 (with a step Ap* = 0.1) and 
several temperatures in the supercritical region. To this 
end, a system consisting of 512 particles was considered. 
The particles were initially placed in a regular configura- 
tion in a cubic volume with periodic boundary conditions, 
with fixed temperature and density. After equilibration, 
the r.d.f. was determined from measurements performed 
over 5 x 10 4 cycles, each of them consisting of an at- 
tempted move per particle. Results for the contact values 
g(l + ) of the r.d.f., as well as for their values g(X~) and 
,g(A+) at both sides of the potential range, were obtained 
from extrapolation and are reported in Table [I] 59 From 
these values, the compressibility factor Z = pV/NkT can 
be obtained from the virial theorem for the SW fluid as 

Z = 1 + l-Kp* {g(l+) - A 3 [g(X-) g(X+)] } . (6) 

Values of Z thus obtained were reported elsewhere^ 8, ex- 
cept for the range A = 1.05. 

TABLE I: MC simulation data of g(l + ), g(A _ ), and g(A + ). 



p' 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 
A = 1.05 

T* = 0.5 

g(l + ) 7.304 7.323 7.343 7.501 7.662 7.771 8.081 8.524 9.078 

g(A~) 7.278 7.250 7.210 7.317 7.374 7.400 7.587 7.832 8.154 

g(A+) 0.987 0.980 0.979 0.988 0.999 0.999 1.030 1.059 1.105 

T* = 0.7 

a(l + ) 4.382 4.609 4.899 5.194 5.557 6.009 6.505 7.223 8.065 

g(A~) 4.321 4.562 4.765 4.986 5.237 5.539 5.842 6.263 6.657 

g(A+) 1.034 1.081 1.142 1.195 1.257 1.325 1.401 1.504 1.595 

T* = 1.0 

g(l+) 2.964 3.215 3.552 3.891 4.336 4.833 5.457 6.258 7.245 

g(A~) 2.923 3.153 3.417 3.679 4.009 4.356 4.759 5.197 5.666 

g(A+) 1.071 1.162 1.255 1.355 1.477 1.599 1.746 1.912 2.085 

A = 1.10 

T* = 0.5 

g(l + ) 7.254 7.080 6.628 6.461 6.088 5.836 5.667 5.620 5.810 

g(A~) 7.179 6.919 6.463 6.219 5.816 5.519 5.307 5.180 5.162 

g(A+) 0.973 0.931 0.875 0.840 0.787 0.746 0.718 0.701 0.698 

T* = 0.7 

g(l+) 4.135 4.153 4.174 4.222 4.334 4.529 4.746 5.064 5.604 

g(A~) 4.097 4.064 4.028 4.007 4.028 4.096 4.135 4.202 4.284 

g(A+) 0.979 0.970 0.964 0.961 0.966 0.980 0.988 1.007 1.027 

T* = 1.0 
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TABLE I: Continued. 



p" 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 

g(l + ) 2.829 2.963 3.134 3.320 3.565 3.868 4.263 4.778 5.481 

g(A _ ) 2.790 2.862 2.974 3.068 3.186 3.316 3.463 3.589 3.713 

g(A+) 1.026 1.058 1.090 1.128 1.172 1.220 1.273 1.316 1.362 

A = 1.20 

T* = 0.7 

g(l + ) 4.126 4.007 3.805 3.562 3.449 3.297 3.281 3.457 3.921 

g(A~) 4.030 3.822 3.555 3.299 3.122 2.924 2.810 2.738 2.612 

g(A+) 0.968 0.914 0.850 0.792 0.748 0.700 0.675 0.657 0.628 

T* = 1.0 

g(l + ) 2.671 2.662 2.666 2.721 2.799 2.941 3.173 3.564 4.232 

g(A~) 2.628 2.551 2.489 2.447 2.415 2.388 2.372 2.326 2.208 

g(A+) 0.968 0.942 0.918 0.900 0.888 0.880 0.871 0.856 0.812 

T* = 1.5 

g(T F ) 2.016 2.092 2.212 2.355 2.541 2.807 3.168 3.692 4.513 

g(A~) 1.956 1.973 1.991 2.019 2.042 2.059 2.063 2.028 1.917 

g(A+) 1.004 1.012 1.021 1.036 1.050 1.057 1.059 1.039 0.982 

A = 1.30 

T* = 1.0 

g(l i ") 2.660 2.592 2.485 2.435 2.414 2.492 2.703 3.136 4.051 

g(A~) 2.584 2.433 2.259 2.138 2.033 1.948 1.865 1.737 1.487 

g(A+) 0.952 0.896 0.830 0.786 0.748 0.716 0.685 0.639 0.547 

T* = 1.5 

g(l + ) 1.957 1.982 2.033 2.105 2.245 2.476 2.829 3.409 4.390 

g(A~) 1.872 1.834 1.792 1.754 1.720 1.680 1.611 1.486 1.272 

g(A+) 0.966 0.940 0.919 0.901 0.883 0.863 0.828 0.763 0.652 

T" = 2.0 

g(l + ) 1.711 1.770 1.880 2.012 2.214 2.497 2.916 3.558 4.537 

g(A~) 1.622 1.618 1.611 1.606 1.582 1.555 1.489 1.368 1.176 

g(A+) 0.990 0.981 0.979 0.972 0.962 0.944 0.903 0.829 0.710 

A = 1.40 

T* = 1.0 

g(l + ) 2.899 2.840 2.640 2.398 2.246 2.245 2.515 3.162 4.351 

g(A~) 2.679 2.488 2.227 1.966 1.783 1.669 1.544 1.332 1.055 

g(A+) 0.986 0.918 0.819 0.725 0.657 0.616 0.568 0.491 0.389 

T* = 1.5 

g(l + ) 1.928 1.916 1.917 1.977 2.075 2.307 2.716 3.442 4.571 

g(A~) 1.836 1.742 1.656 1.583 1.518 1.441 1.328 1.149 0.932 

g(A+) 0.945 0.894 0.849 0.811 0.778 0.740 0.681 0.589 0.478 

T* = 2.0 

g(l + ) 1.665 1.705 1.788 1.889 2.070 2.362 2.834 3.584 4.719 

g(A~) 1.588 1.533 1.489 1.444 1.400 1.331 1.227 1.067 0.866 

g(A+) 0.960 0.929 0.905 0.877 0.852 0.807 0.743 0.646 0.525 

A = 1.50 

T* = 1.5 

g(l + ) 1.952 1.909 1.884 1.888 1.989 2.263 2.783 3.640 4.962 

g(A") 1.832 1.695 1.570 1.464 1.382 1.281 1.150 0.993 0.850 

g(A+) 0.941 0.870 0.804 0.751 0.709 0.659 0.590 0.510 0.437 

T* = 2.0 

g(l + ) 1.661 1.660 1.720 1.820 2.006 2.336 2.880 3.740 5.013 

g(A~) 1.563 1.478 1.403 1.338 1.274 1.185 1.063 0.922 0.790 

g(A+) 0.945 0.894 0.851 0.814 0.771 0.719 0.646 0.558 0.480 

T* = 3.0 

g(l + ) 1.444 1.521 1.628 1.795 2.046 2.424 2.998 3.825 5.025 

g(A~) 1.351 1.313 1.268 1.230 1.171 1.090 0.982 0.857 0.728 

g(A+) 0.968 0.938 0.909 0.878 0.839 0.781 0.704 0.614 0.521 

A = 1.60 

T* = 1.5 

g(l + ) 2.064 2.026 1.949 1.894 2.004 2.339 3.002 3.996 5.382 

g(A~) 1.869 1.696 1.528 1.387 1.293 1.200 1.099 1.015 0.956 

g(A+) 0.958 0.872 0.784 0.713 0.665 0.618 0.565 0.521 0.491 

T* = 2.0 

g(l + ) 1.655 1.673 1.702 1.795 2.007 2.393 3.037 3.985 5.313 



TABLE I: Continued. 



p* 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 

g(X~) 1.542 1.446 1.350 1.271 1.199 1.110 1.015 0.932 0.872 

g(A+) 0.940 0.873 0.819 0.770 0.726 0.673 0.615 0.564 0.530 

T* = 3.0 

g(l + ) 1.439 1.498 1.602 1.762 2.040 2.463 3.087 3.987 5.258 

g(A~) 1.329 1.269 1.218 1.165 1.100 1.021 0.933 0.850 0.792 

g(A+) 0.953 0.912 0.873 0.834 0.788 0.732 0.669 0.611 0.567 

A = 1.70 

T* = 2.0 

g(l^) 1.716 1.714 1.752 1.833 2.077 2.538 3.250 4.253 5.611 

g(A~) 1.551 1.432 1.322 1.237 1.165 1.101 1.045 1.017 1.022 

g(A+) 0.944 0.870 0.802 0.749 0.706 0.668 0.634 0.617 0.619 

T* = 3.0 

g(l + ) 1.430 1.494 1.598 1.779 2.083 2.549 3.222 4.160 5.454 

g(A~) 1.318 1.251 1.184 1.131 1.069 1.009 0.956 0.928 0.933 

g(A+) 0.949 0.893 0.849 0.810 0.766 0.725 0.685 0.664 0.667 

T* = 5.0 

g(l + ) 1.299 1.403 1.554 1.779 2.114 2.567 3.220 4.108 5.351 

g(A~) 1.177 1.139 1.097 1.052 1.000 0.941 0.893 0.860 0.867 

g(A+) 0.964 0.934 0.898 0.861 0.818 0.771 0.731 0.703 0.709 

A = 1.80 

T* = 2.0 

g(l + ) 1.839 1.916 1.880 1.932 2.199 2.708 3.487 4.508 5.886 

g(A~) 1.599 1.468 1.326 1.233 1.175 1.144 1.141 1.183 1.312 

g(A+) 0.974 0.892 0.805 0.745 0.714 0.694 0.693 0.716 0.793 

T* = 3.0 

g(l + ) 1.449 1.522 1.629 1.828 2.162 2.670 3.392 4.325 5.627 

g(A~) 1.320 1.238 1.175 1.124 1.079 1.046 1.043 1.079 1.195 

g(A+) 0.942 0.888 0.840 0.803 0.773 0.751 0.747 0.772 0.853 

T* = 5.0 

g(l + ) 1.297 1.411 1.562 1.809 2.152 2.646 3.318 4.214 5.465 

g(A~) 1.169 1.126 1.085 1.042 1.005 0.978 0.970 1.003 1.107 

g(A+) 0.956 0.921 0.886 0.854 0.822 0.800 0.794 0.821 0.906 

A = 1.90 

T* = 3.0 

g(l + ) 1.492 1.570 1.692 1.908 2.274 2.815 3.533 4.450 5.534 

g(A~) 1.323 1.247 1.182 1.142 1.123 1.132 1.174 1.265 1.403 

g(A+) 0.949 0.895 0.848 0.818 0.805 0.811 0.841 0.908 1.008 

T* = 5.0 

g(l + ) 1.309 1.416 1.587 1.847 2.222 2.733 3.404 4.268 5.353 

g(A~) 1.170 1.123 1.087 1.058 1.044 1.049 1.089 1.174 1.305 

g(A+) 0.958 0.920 0.889 0.867 0.855 0.859 0.891 0.961 1.072 

A = 2.00 

T* = 3.0 

g(l+) 1.564 1.699 1.824 2.032 2.424 2.973 3.669 4.493 5.499 

g(A~) 1.346 1.282 1.217 1.192 1.199 1.245 1.322 1.425 1.528 

g(A+) 0.968 0.918 0.873 0.853 0.861 0.892 0.949 1.024 1.100 

T* = 5.0 

g(l + ) 1.325 1.450 1.633 1.907 2.300 2.814 3.471 4.304 5.354 

g(A~) 1.178 1.131 1.105 1.098 1.109 1.152 1.225 1.326 1.428 

g(A+) 0.960 0.926 0.905 0.899 0.910 0.944 1.005 1.089 1.175 



IV. RESULTS AND DISCUSSION 

The objective is to determine the limits of applicabil- 
ity of the TL and YS theories. We are mainly inter- 
ested in the domain of moderate temperatures. By that 
we mean temperatures within the range T*(X) £ T* ^ 
3T*(A), where T*(X) denotes the critical temperature of 
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FIG. 1: The open circles represent the values of the reduced 
temperature T* we have considered for each value of the range 
A. The crosses are simulation data for the critical tempera- 
ture T c * (A) (see Ref. I24I) and the solid line is the theoretical 
estimate Jjj. Note the logarithmic scale of the vertical axis. 



the SW fluid with range A. This critical temperature 
has been measured in computer simulations for several 
rangesi 39 ! 40 ' 41 ! 42 ! 43 ! 44 ! 53 Table I of Ref. |H gives a rather 
extensive compilation of data. A simple analytical esti- 
mate for T*(A) was derived in Ref. l30t 
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FIG. 2: Comparison of the contact values g{l + ) of the r.d.f. 
obtained from the TL (dashed lines) and YS (solid lines) theo- 
ries with Monte Carlo data as functions of the reduced density 
p* for different temperatures and well widths A = 1 .05-1 .5. 



T*(A) = 



In 



i , 3+A+2V2A 
1 A(A-1)(9-2A+A a ) 



(7) 



FigureGJis a T*-A plot where the open circles represent 
the three temperatures (two in the cases A = 1.9 and 
A = 2) we have considered in the simulations for each 
value of A. The simulation data of T*(A)^4 as well as the 
theoretical estimate are also shown. For each value of 
A we have typically considered three temperatures: T* > 
T*, T* « 1.5T*, and T* « 2T*, so that T* > 3T C *(A). 

Before comparing the simulation data of the full r.d.f. 
with the theoretical predictions, it is worth focusing on 
the contact value fl'(l + ). Figures and show g(l + ) as 
a function of the reduced density for the temperatures 
represented in Fig. ^ as obtained from our MC simula- 
tions, as well as from the TL and YS theories [cf. Eqs. 
lUD, (SHI, and {B5}]. We observe that for A = 1.05 and 
A = 1.1 the non-perturbative YS model presents a very 
good agreement with the simulation data for the three 
temperatures considered, whereas the TL perturbation 
theory is rather poor, especially for low temperatures. 
For A > 1.2, however, the YS model behaves well for 
small and moderate densities but starts to fail in the 
high-density domain, especially for the lowest tempera- 
ture, the failure being more dramatic as the well width 
increases. Interestingly, the TL theory becomes more 
accurate precisely in that high-density region where the 
YS model is less reliable. Thus, for a given range A, 
there exists a certain threshold density pg(X) such that 
the YS model is accurate for p* < Po(A) an d inaccurate 



for p* > Pq(A), while the opposite situation occurs in 
the case of the TL theory. Of course, this qualitative 
description applies for the range of "moderate" tempera- 
tures defined above, since the results obtained from both 
theories tend to coincide as the temperature increases. 

According to Figs. |3 and |3| the location of Pq(X) 
roughly coincides with the region where either the 
isotherms cross (for A < 1.7) or have the least separation 
(for 1.8 < A < 2.0). This means that the simulation data 
of ff(l + ) in the region p* w /?q(A) are practically insen- 
sitive to the temperature, so they are close to its hard- 
sphere value 5o(l + )- For larger densities, p* > Pq(X), 
the simulation data show that the influence of tempera- 
ture is small and hence the TL perturbation theory be- 
comes accurate in that domain. On the other hand, for 
P* PoW the MC values of g(l + ) are strongly sensi- 
tive to temperature, as expected from the fact that at 
zero density <?(1 + ) = e 4 / T , while perturbation theories 
give g(l + ) = l + l/T*. The strong deviation of the non- 
perturbative YS theory from the MC data in the density 
region p* > Po(A), especially for A > 1.4, is in part due 
to the fact that in the YS model the parameter A in 
Eq. Q is assumed to be independent of density and so 
it is assigned its zero-density value A — e 4 / T — 1. A 
better agreement is expected if A is allowed to depend 
on density, but this would imply either to impose an ex- 
tra consistency condition (for instance, continuity of the 
first derivative of the cavity function) or to apply an em- 
pirical fit, what is outside the original spirit of the YS 
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FIG. 3: Same as in Fig.H but for A = 1.6-2. 
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FIG. 4: Plot ol the threshold density Po(A) as a function of 
the potential range A. For each value of A, po is defined as 
the density around which the MC contact values g(l + ) are 
practically insensitive to the temperature. Below (above) the 
curve, the YS (TL) theory can be considered as reliable. The 
line is a guide to the eye. 
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FIG. 5: Comparison of the r.d.f. obtained from the TL (dot- 
ted lines) and YS (solid lines) theories with Monte Carlo data 
(circles) for A = 1.05 and T* = 0.5. Note that the TL curves 
are interrupted after x = 2. 



model. A second reason has to do with the construction 
of the YS model as an extension of Baxter's solution of 
the PY integral equation for sticky hard spheres, so that 
in principle it is intended to be a model for narrow wells. 

A plot of Pq(X) is presented in Fig. It can be inter- 
preted as a sort of "phase" diagram in which the curve 
separates the respective regions where the TL and YS 
theories are reliable for moderate temperatures in the in- 
terval T*(A) < T* < 3T*(A). We observe that as the 
range A decreases, the YS region tends to span the whole 



fluid density domain. In addition, pg presents a mini- 
mum Pq w 0.4 at A « 1.7, so the YS theory does a fairly 
good job if p* < 0.4, even for wide potentials. 

Once we have analyzed the performances of the TL and 
YS theories in connection with the contact value 
let us proceed to investigate the r.d.f. g(x) itself. The 
results are presented in Figs.lS HTUI Since in this paper we 
are mainly interested in short-ranged SW potentials, we 
have paid special attention to the ranges 1.05 < A < 1.3 
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FIG. 6: Same as in Fig. 03 but for A = 1.1 and T* = 0.5. 
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FIG. 7: Same as in Fig. ED but for A = 1.2 and T* = 0.7. 



(Figs. EDS. As representative examples of a moderate 
and of a wide range we have considered A = 1.5 (Fig. 
El) and A = 2.0 (Fig. 110(1 . respectively. For each value of 
A we have restricted ourselves to the lowest temperature 
represented in Fig. ^ and to the densities p* — 0.2, 0.4, 
and 0.8 (except in the case A = 2.0, where p* = 0.8 has 
not been considered because the YS model fails to have a 
solution in that case). In agreement with the analysis of 
Figs. [21 and E3 one can see that the YS theory works well 
for small potential widths (A < 1.2) for the whole density 
range. For larger potential widths, the performance of 
the theory is still fair at low (p* = 0.2) and even moderate 



(p* = 0.4) densities. However, the YS theory fails, and 
even can become entirely unphysical, at high densities 
(p* = 0.8 > pg). Of course, at temperatures higher than 
those of Figs. IKlllOlthe performance of the theory at high 
densities improves (not shown). 

By contrast, the TL theory presents the opposite be- 
havior to that of the YS theory since its accuracy in- 
creases as the density and the potential width grow. Of 
course, it also improves if the temperature increases, as 
expected from a perturbation theory. According to Figs. 
IS MlQI the TL theory does a better job than the YS model 
at p* — 0.8 for A > 1.3, in agreement with the "phase" 
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diagram of Fig. 0] 

V. CONCLUSIONS 

In this paper we have presented extensive Monte Carlo 
simulations for the structural properties of square-well 
fluids with ranges A, reduced densities p* , and reduced 
temperatures T* in the intervals 1.05 < A < 2, 0.1 < 
p* < 0.8 and T*(A) < T* < 3T*(A), respectively. The 
MC data have been used to assess the accuracy of two 
theories that provide explicit expressions of the r.d.f. in 
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Laplace space, the TL perturbation theory2L22i and the 
non-perturbative YS models 

The results show that both theories complement each 
other, as the YS theory works well where the TL theory 
fails and vice versa. More specifically, the YS theory ex- 
hibits a good agreement with the MC data at any fluid 
density if the potential well is sufficiently narrow (say 
A < 1.2), as well as for any width if the density is small 
enough (say p* < 0.4). This can be further refined by 
noticing that the YS theory works well if p* < Po(^)' 
where Pq(X) is the density around which the simulation 
data for the contact value g(l + ) show the least influence 
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APPENDIX A: EXPLICIT EXPRESSIONS IN 
THE TANG-LU THEORY 

Let us define the Laplace transform G(s) of xg(x): 

/oo 
dxe~ tx xg{x). (Al) 

The contact value <?(1 + ) is given from G(s) as 

, 9 (1+) = lim t^Gtt). (A2) 

t — >oo 

The exact solution of the PY equation for hard 
spheres^*^ reads 



G (t)=t 



L(t)e- 



S(t) + 12r)L(t)e- 



where i] — §pcr 3 is the packing fraction and 



L(t) = 1 + 2i] + (1 + tj/2)t, 



(A3) 



(A4) 



on temperature. On the other hand, for p* > Pq(X) the 
YS theory rapidly deteriorates, especially for tempera- 
tures near the critical one, while the TL theory becomes 
very accurate. 



The complementarity between the TL and YS theories 
is interesting because they present some formal similari- 
ties in their formulation and are (practically) equally easy 
to implement (see Appendices lAl and |B"|) . The latter the- 
ory, however, has some advantages over the former one. 
First, the YS theory is especially useful for describing 
colloidal dispersions modeled as short-ranged SW fluids. 
Second, it provides a simple analytical expression for the 
second shell (2 < x < 3) of the r.d.f., whereas this is not 
the case for the TL theory. 28 Last, it seems feasible to im- 
prove the performance of the YS theory at high densities 
by imposing additional constraints to the Laplace trans- 
form of the r.d.f. to determine the parameter A in Eq. 
as a function of density. Instead, in order to improve 
the TL theory it would be necessary to obtain higher or- 
der terms in the expansion of the r.d.f. of the SW fluid 
in power series of the inverse of the reduced temperature 
T*, and this seems to be too complicated at present. 



S{t) = -12^(1 + 2rj) + 18r7 2 t + 677(1 - V )t 2 + (1 - r]) 2 t 3 . 

(A5) 

The corresponding contact value is 



3o(l + ) = 



1 + v/2 

(1-77)2- 



(A6) 



Equation l|A3(l provides the zeroth-order term in the 
TL perturbation theory. The first-order term is2£*2& 



Gi(t) 



(l-„)«e-* r * 4 (l + At) ^ A _^ 



?§(*) 



s 2 (-t) 



, S 2 (U) 
'Si{U) 



t + t. 

4 + (1 + 4A)tj + A (A - l)tf 



+ ti(l-Xti) 



xe 



(A- 



-1)**} 



(A7) 



where Si(t) = S'(t), S 2 (t) = S"(t), the primes denoting 
derivatives with respect to t, and 



Qo(t) 



S(t) + l2r)L(t)e 

(1 -T]) 2 t 3 



(A8) 
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In Eq. i|A7|l the summation extends over the three zeroes 
of S(t), denoted by fcj. The contact value gi(l + ) is 



91 



=J sf(u) 



Si(U) 



+(1 + 4A)tj + A(A - l)tf\ e^ 1 ^ 



(A9) 



By analytical inversion of G\{t) one can get explicit 
expressions for gi(x) inside the shells n < x < n + 1, 
which become increasingly more complicated as n grows. 
The expression for the first shell 1 < x < 2 can be found 
in Ref.ES 



Here, ij are the three distinct roots of S(t) = 1 + Sit + 
S 2 t 2 + S3t 3 [not to be confused with the polynomial i|A5[l ] 
and 

C w (t) = l + A + Kit, C n (t) = -(A + K 2 t). (B8) 
Analogously, the inverse Laplace transform of t[F(t)] 2 

is 

f 2 (x) = f w {x)Q{x)+f 2 i(x + l-\)Q(x + l-\) 



-f 22 (x + 2-2X)e(x + 2-2X), 



where 



(B9) 



APPENDIX B: EXPLICIT EXPRESSIONS IN 
THE YUSTE-SANTOS MODEL 

By imposing the exact condition G(t) — t~ 2 ~ t for 
small t, where G(t) is defined by Eq. IjAljl . one can ex- 
press the parameters K%, Si, S 2 , and S3 appearing in Eq. 
(QJ as linear functions of A and if 2 i 29 i 30 



1 + 2?/ 



1 + 1 + 2r/(A 3 - l)K 2 



-^(A 4 - 4A + 3)Aj + K 2 - A(X - 1), (Bl) 



Si 



l + 2?7 
1 



| + 2(A 3 -l)if 2 -i(A 4 -4A + 3)^ 



{—1 + 77 + 2 [A — 1 — 2??A(A 2 - 1)] K 2 

(B3) 



2(1 + 2?]) 

-[(A-l) 2 -rKA 2 -l) 2 ]A}, 



S3 = 



1 



-(A 2 -l)- ?/ A 2 (A-l) 



K 2 



c (t) s " iu) 



a tiX 



where we have called 

C 20 (t) = t[C w {t)] 2 , C 21 (t) = 2tC w (t)C n (t), 
C 22 (t) = t[Cn(t)} 2 . 



(BIO) 



(Bll) 



Insertion of Eqs. (|B6f) and l|B9|l into Eq. JSJ gives the 
r.d.f. g(x) in the interval 1 < x < 3. Note that the 
contribution f 22 {x) is needed inside that interval only if 
A < |. For x > 3 the evaluation of f 3 {x), f±(x), ... is 
required. Alternatively, one can make use of the effi- 



(B2) cient method discussed by Abate and WhitlpSi to invert 



(I-??) 2 
1 + 2ri { 12rj 

+ ^[4 + 2A - r;(3A 2 + 2A + 1)](A - I) 2 A ) . (B4) 



numerically Laplace transforms. 

To close the model, we need to determine the param- 
eters A and K 2 . The former is assigned its zero-density 
limit value, namely A = e x l T — 1»S To determine K 2 we 
impose the continuity condition of the cavity function at 
x = A, what implies 



<?(A- 



This yields 



From Eq. I|A2(1 . we have 



9(l + ) = 



_Ki_ 

12 V S 3 



(l-e- 1/T *)/ 10 (A-l) = -/ 11 (0) = - 



12 V S 3 



(B12) 



(B13) 



(B5) 



By application of the Heaviside expansion theorem, the 
inverse Laplace transform of tF(t) reads 

fx{x) = fio(x)e(x) + f n (x + 1 - A)9(x + 1 - A), (B6) 

where 

3 



r 1 x 1 \ " Cikjtj) + 

Ilh[X) ~ 12 V ^ S'(ti) 1 



(B7) 



Since the roots U depend on K 2 through the coefficients 
Si, S 2 , and S3, Eq. I|B13|I is a transcendent equation for 
K 2 that needs to be solved numerically. Acedo and San- 
tos have recently proposed a simplified version of the YS 
model whereby the exact condition i|B12(l is replaced by a 
simpler one that allows K 2 to be obtained analytically. 30 
This is especially useful for determining the thermody- 
namic properties In this paper, however, since we 
are interested in the structural properties, we enforce 
condition (|B12|) and determine K 2 from Eq. (|B13I) . 
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